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Abstract 

The skyrmion is a topologically stable spin texture, in which the spin direction wraps 
a spherei. The topological nature gives rise to emergent electromagnetic phenomena 
such as topological Hall effect^"^. Recently, the crystallization of nanoscale skyrmions 
was observed in chiral helimagnets with use of the neutron diffraction and Lorentz 
transmission electron microscopy^ The skyrmions are quite mobile under electric 
current density as low as ~ 10^ A/m^, and, in some cases, stable up to near room 
temperature^'^. These features suggest that the skyrmions may work as a magnetic 
stable variable similar to the bubble memory— but equipped with topological func- 
tionalities. Here, we investigate the low-energy excitations of the skyrmion crystal 
in a helimagnetic insulator Cu20Se03^° in terms of microwave response. We have 
observed two distinct excitations of the skyrmion with different polarization charac- 
teristics; the counter-clockwise circulating mode at 1 GHz and the breathing mode 
at 1.5 GHz. These modes may play a crucial role in the low energy dynamics of 
skyrmions and hence in their manipulation via external stimuli. 
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The manipulation of nano- or micro-scale magnetic objects has been the central issue 
of spin-related electronic science and technology. The bubble domains in the ferromagnets 
with uniaxial magnetic anisotropy can easily be manipulated by magnetic field gradient, and 
therefore were applied to the memory storage devise^. The electrical control of ferromagnetic 
domain wall is thought to be a key technology for the magnetic random access memoryii. 
Recently, a topological magnetic object termed skyrmion has been attracting much attention 
since the triangular lattice of skyrmions was observed in chiral helimagnets such as MnSi and 
(Fe,Co)Si^i^. As shown in Fig. Ic, the magnetic moments composing the skyrmion point 
downward (opposite to the magnetic field) at the core and upward on the edge, and are 
aligned circularly in between. It may be viewed as a nanoscale version of bubble memory 
but has a unique topological nature. Because the direction of magnetic moment varies 
continuously and wraps all the solid angle in the skyrmion, the skyrmion carries a topological 
quantum number termed skyrmion number 



where n is the unit vector parallel to the magnetic moment at the position r = {x,y,z). 
Reflecting the topological nature, the skyrmion works as a fictitious magnetic flux quantum 
4>o = h/e and gives rise to emergent electromagnetic phenomena such as topological Hall 
effect, spin- motive force, and skyrmion Hall effect^~-»^>^. 

The important issue to be clarified is the dynamics of skyrmions. It has recently been 
demonstrated that the skyrmion can be easily driven by electric current^iiSi. Jonietz et 
al. reported that the rotation of diffraction pattern is caused by the electrical current in 
the presence of the thermal gradient and ascribed it to the result of the current drive of 
skyrmions^. On the other hand, Schultz et al. have observed the emergent electric field or 
spin motive force induced by the skyrmion current^. The threshold electric current (~ 10^ 
A/m^) is five orders of magnitude lower than that needed to drive the ferromagnetic domain 
wall, which implies the potential of practical application. 

The resonant excitation by microwave irradiation provides a useful probe for dynamics 
of magnetic objects. In fact, the magnetic bubbles can be created and detected with use 
of the microwave^. The oscillation of magnetic domain walls can also be excited by the 
AC magnetic field with the frequency of sub-GHz ranged. Here, we have investigated the 
magnetic excitations of skyrmion crystal with use of microwave in a helimagnetic insulator 
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Cu20Se03 with small Gilbert damping. 




FIG. 1: Magnetic excitations in skyrmion crystal and helimagnetic state, a, Magnetic 
phase diagram for Cu20Se03 as determined by the magnetic susceptibiUty, which is measured 
with minimal demagnetization field (corresponding to the Hbc\\Hac configuration in the following 
microwave experiments), b, Real-space image of a skyrmion in Cu20Se03 obtained by Lorentz 
transmission electron microscopy, reproduced from ref. 10. c-f, Illustrations of skyrmion (c) and 
magnetic excitations in skyrmion crystal (d-f). g., Illustrations of magnetic excitations in helim- 
agentic state. The untransparent and half transparent magnetic moments indicate the magnetic 
moments in static and excited states, respectively. Bottom panel illustrates the oscillating compo- 
nent (the difference between the excited and static states) of magnetic moments in the coordinate 
system fixed to the static magnetic moment. The numbers stand for the positions of magnetic 
moments. 

While the skyrmion crystal has been observed so far mostly in the B20 transition metal 
compounds such as MnSi and (Fe,Co)Si, it has recently been found that it can be realized 
also in an insulating oxide Cu20Se03^*^. In common with the B20 transition metal com- 
pounds, the space group of the crystal structure is P2i3, which is noncentrosymmetric but 
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nonpolar and cubic. Figure la shows the magnetic phase diagram for Cu20SeO^. The 
S = 1/2 moments at the Cu^"*" sites show the local ferrimagentic arrangement of three- 
up and one-down type below T^, ~ 60 K as observed by power neutron diffraction^^ and 
NMR— measurements. However, the Dzyaloshinskii-Moriya interaction arising from the 
noncentrosymmetric lattice structure modulates the ferrimagnetic correlation and induces 
the helical spin structure with the period of about 50 nm at zero magnetic field. While the 
magnitude of the helical wave vector is determined by the ratio between the asymmetric 
Dzyaloshinskii-Moriya and symmetric ferrimagentic interactions, the wave vector direction 
is nearly degenerate, and the multi-domain structure is formed at zero magnetic field^S'. 
When the magnetic field is applied, the wave vector tends to be aligned along the field di- 
rection, forming the single-domain conical spin structure. In a high enough magnetic field, 
the induced ferrimagnetic state is realized. Just below the critical temperature, on the other 
hand, there appears the nontrivial magnetic state, i.e. skyrmion crystaU^. The real space 
image of skyrmion as obtained by Lorentz transmission electron microscopy is reproduced 
in Fig. lb, where the in-plane components of magnetic moments in the skyrmion are well 
resolved. 

The skyrmion crystal in a bulk crystal is restricted to a narrow range of magnetic field, 
therefore the microwave response can hardly be observed in a conventional ESR spectrome- 
ter, in which the probe frequency is fixed and the magnetic resonance is probed by scanning 
magnetic field. Thus, we have constructed a broad-band microwave system with a network 
analyzer as detailed in the Supplementary Information. While the pioneering work of elec- 
tron spin resonance for MnSi was previously done with use of several fixed frequencies^^, 
our broad band measurement could provide much more information as shown below. At 
first, we show the microwave response at 40 K to discuss the magnetic excitations in the 
helimagnetic state (Scan 1 shown in Fig. la). In Figs. 2a and b, we show the magnetic 
field dependence of microwave absorption spectra A^n at 40 K. (Here, ASn stands for 
the difference of the reflection coefficient 6*11 of the sample-inserted transmission line from 
the background value, which is proportional to the microwave absorption of the sample; for 
details, see Supplementary Information). The AC magnetic field of microwave i^Ac is per- 
pendicular to the DC magnetic field -ffoc in a, and parallel in b. For the both cases, Hy)q is 
applied along < 110 > direction. The microwave spectra are qualitatively unchanged by the 
variation of i^DC direction with respect to the crystal axis although they are affected more 
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FIG. 2: Magnetic field dependence of microwave absorption spectra at 40 K. a,b, Mi- 
crowave absorption spectra A5ii at various magnetic fields for H^c ± ^ac (a) and Hdc\\Hac (t»). 
c, The magnetic field variation of magnetic susceptibility at 40 K measured in the configuration 
corresponding to the Hdq -L Hac microwave measurement, d The magnetic field dependence of 
the resonance peak frequency in the microwave wave absorption spectra for Hdq _L -f^AC- 



or less by the static demagnetization field since we used a plane-like sample in order to fit 
it to the microwave probe. In the configuration of i^DcH-f^AC! -f^DC is parallel to the sample 
plane and the demagnetization field is small. On the other hand, -ffoc is perpendicular to 
the sample plane and the demagnetization field is large in the ifnc -L -f^Ac configuration. 
A^ii shows broad two peaks around 4 GHz for both ifod |-f^Ac and i^DC -L -f^Ac in the low 
i^DC region. The peaks increase in intensity and become sharper when i^DC = 400 Oe is 
applied perpendicular to H^c whereas they are completely suppressed above 200 Oe in the 
parallel configuration. As shown in Fig. la, the magnetic wave vector q is aligned along 
Hy)C forming a single g-domain above 200 Oe. (For the H^c -L -f^Ac configuration, the 
critical magnetic field is ~ 400 Oe becasue of the large demagnetization field.) Therefore, 
the contrastive polarization dependence in the high magnetic field indicates that the two 
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magnetic modes can be excited only by Hac perpendicular to the helimagnetic wave vector. 
Kataoka theoretically identified such two modes in the helical magnetic state induced by 
Dzyaloshinskii-Moriya interaction^^. The two modes, +Q mode and -Q mode, are shown in 
Fig. Ig. In the coordinate fixed to the static magnetic moment direction, the oscillating 
component of +Q mode at a fixed time rotates along the helical wave vector similarly to the 
static helical moments as illustrated in the bottom left panel of Fig. Ig. On the other hand, 
in the -Q mode, the rotating direction is opposite. The +Q and -Q modes are almost degen- 
erate but the degeneracy may be lifted by the additional antiferromagnetic interactioniS. As 
for the magnetic field dependence in the high field region above 400 Oe, the frequencies of 
observed two peaks decrease with i^oc below 1000 Oe. Above 1500 Oe, only one peak is ob- 
served, whose frequency rather increases with i^DC- Iii Figs. 2 c and d, we show the magnetic 
susceptibility x and the frequency of the observed peaks in the perpendicular configuration 
as functions of magnetic field. The magnetic susceptibility is measured in magnetic fields 
perpendicular to the sample plane corresponding to the i^DC -L Hac configuration. The 
metamagnetic transition from the conical magnetic state to the induced ferrimagnetic state 
is observed around 1400 Oe as a kink of magnetic susceptibility. As indicated by a thick 
gray bar, the change in the observed peak magnetic resonance mode(s) and its(their) field- 
dependence corresponds to this metamagnetic transition. The magnetic -field dependence of 
the mode frequency is in accord with the theoretical resullr^, thereby confirming the known 
character of the magnetic excitations in spin helix. 

Next, we show the newly resolved magnetic resonances in the skyrmion. We have inves- 
tigated the magnetic field dependence of microwave absorption spectra at 57.5 K to probe 
the magnetic excitations in skyrmion crystal (Scan 2 shown in Fig. la). In Fig. 3a, we 
plot AS'ii for i^DC -L -f^Ac at various magnetic fields. In the low-field helical magnetic state, 
one broad excitation is observed around 1.8 GHz. The splitting is not clearly observed in 
this temperature and the frequency is lower than that at 40 K. The peak is unchanged or 
slightly enhanced with i^oc up to 100 Oe. Nevertheless, the helical peak becomes quite 
weak and alternatively another peak emerges around 1 GHz above 140 Oe. Around 340 Oe, 
the low frequency peak disappears and the helical mode revives. In Figs. 3b-d, we show 
the frequency and intensity of the peaks as functions of magnetic field, comparing them 
with the magnetic susceptibility measured in this configuration. The susceptibility shows 
peaks around 120 Oe and 350 Oe, indicating the phase boundary between the helimagnetic 
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FIG. 3: Microwave response for i?AC -L ^^DC in skyrmion crystal, a, The microwave 
absorption spectra A5ii for Hdq -L -f^AC at various magnetic fields at 57.5 K. b, The magnetic 
field variation of magnetic susceptibility at 57.5 K measured in the configuration corresponding to 
the Hdc _L -ffAC microwave measurement. c,d, The frequency (c) and intensity (d) of the magnetic 
modes in the microwave absorption spectra in a. e, The intensity of the magnetic-resonance mode 
for /^DC -L ^AC (counter clock wise rotaional mode) in skyrmion crystal plotted in the H — T phase 
diagram determined by the magnetic susceptibility measurement. Squares stand for the measured 
(temperatures, magnetic field) points. 

and skyrmion crystal states. The field range, where the low-frequency peak is observed, 
coincides with that for the skyrmion crystal phase. The intensity of the helical mode is 
not completely suppressed even in the skyrmion crystal phase. A similar coexistence of the 
helical and skyrmion crystal states is hinted by the neutron diffraction measurement for B20 
compounds (the coexistence of the diffraction spot of q\\H due to the helical state and that 
of g _L if due to the skyrmion crystal)^. In the induced ferrimagnetic state above 600 Oe, 
the magnetic susceptibility is suppressed and the peak frequency readily increases with Hdq. 
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Figure 3e shows the contour mapping of the low-frequency peak intensity in the H — T phase 
diagram around the skyrmion-crystal phase determined by the magnetic susceptibihty mea- 
surement. It is clear from this figure that the low-frequency mode at 1.0 GHz is emergent 
only in the skyrmion crystal phase. Thus, the low-frequency peak is assigned to a magnetic 
excitation inherent in the skyrmion crystal. 
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FIG. 4: Microwave response for Hxq\\Hj)q in skyrmion crystal, a, The magnetic field 
dependence of microwave absorption spectra A^n for Hdq\\Hxq at 57.5 K. b, The magnetic field 
variation of magnetic susceptibility at 57.5 K measured in the configuration corresponding to the 
-f^DcII^AC microwave measurement. c,d, The frequency (c) and intensity (d) of the magnetic mode 
in the microwave absorption spectra for H-£,q\\Hxc- e, The intensity of the magnetic-resonance 
mode for H-qqWH^q (breathing mode) in skyrmion crystal plotted in the H — T phase diagram. 
Squares stand for the measured (temperatures, magnetic field) points. 



Figure 4a shows the magnetic field variation of AS*!! for i^DcH-f^Ac at 57.5 K. The mag- 
netic excitation cannot be excited by the AC magnetic field parallel to the magnetic moment. 
Note here that, conventional magnetic excitations are usually not observed in this configura- 
tion. Nevertheless, we discerned a peak around 1.5 GHz between 50 Oe and 150 Oe. In Fig. 
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4b-d, we plot the frequency and intensity of the peak as well as the magnetic susceptibility 
measured in this configuration. In the magnetic susceptibility, the sharp peaks are observed 
at 75 Oe and 250 Oe indicating the phase transition between the conical magnetic phase 
and skyrmion crystal phase. The reason why the magnetic field region of skyrmion-crystal 
phase is lower than that in the -f^DC -L -f^Ac configuration is because of the demagnetization 
field as mentioned previously. The peak around 1.5 GHz is emergent around this magnetic 
field range. Also for this Hac\\H-dc configuration, we show the contour mapping of the 
peak intensity in the magnetic phase diagram determined by the magnetic susceptibility 
measurement. The region, where the absorption peak is emergent, almost coincides with 
the skyrmion crystal phase. (The small difference of the magnetic field (~ 50 Oe) is caused 
by the artifact perhaps due to the flux pining in the superconducting magnet.) The mi- 
crowave absorption peak at 1.5 GHz in this parallel configuration can also be assigned to 
the excitation in the skyrmion crystal. 

Recently, Mochizuki numerically identified three spin-wave modes in skyrmion crystal 
on the basis of Landau- Lifshitz-Gilbert equation^^. Two modes are counter-clockwise and 
clockwise rotational modes illustrated in Figs. Id and e, in which the core of skyrmion 
rotates in the counter-clockwise and clockwise directions, respectively. These modes can 
be excited by the in-plane Hac- Petrova and Tchernyshyov analytically derived the similar 
rotational modes^. The third mode is the breathing mode shown in Fig. If, in which the core 
of skyrmion expands and shrinks alternately and can be excited by the out-of-plane H^c- 
Since the skyrmion plane is always perpendicular to the DC magnetic field, the experimental 
configurations of Hac -L -f^DC and Hac\\Hbc correspond to the in-plane and out-of-plane 
Hac, respectively. The observed mode at 1.5 GHz in the Hac\\Hbc configuration (Fig. 
4) can be undoubtedly assigned to the breathing mode of the skyrmion. For the in-plane 
polarization {Hac -L H^q), on the other hand, only one skyrmion mode is experimentally 
observed, while two rotational modes are theoretically anticipated. It is, however, to be noted 
that the intensity and frequency of the counter-clockwise rotational mode are larger and 
lower, respectively, than those of the clockwise rotational mode according to the theoretical 
simulation. The experimentally observed mode corresponds perhaps to the intense counter- 
clockwise rotational mode, while the higher lying clockwise mode may be mixed with the 
remaining helical mode and not clearly identified in the experiment. 

In summary, we have observed the breathing and rotational modes of skyrmion in a 
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helimagent Cu20Se03. The internal vibrations of skyrmion should be endowed with the 
magnetoelectric coupling. Therefore, nontrivial microwave responses such as directional 
dichroism are expected similarly to electromagnons in the THz regime^i. In addition, these 
modes will be useful for the manipulation of skyrmions as already demonstrated by the 
numerical calculation^^. 

Methods 

Single crystals of Cu20Se03 were grown by the chemical vapor transporlr^i^. Magnetic 
susceptibility was measured in a SQUID magnetometer. The broad band microwave ab- 
sorption spectra was measured in a superconducting magnet with use of a vector network 
analyzer. The detail of the microwave experiment is given in the Supplementary Information. 
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Supplementary Information:Experimental setup of microwave absorption 

In order to observe magnetic excitations in skyrmion crystal, we have set up a broad band 
microwave measurement system with a vector network analyzer. The experimental setup 
is depicted in Fig. 5a. The semi rigid cable connected to the network analyzer is inserted 
into a superconducting magnet. A microstrip line^ composed of narrow and broad copper 
plates is attached at the end of semi rigid cable (Fig. 5b). The other end of microstrip 
line is electrically shorted. Because the characteristic impedances of the semi rigid cable 
and microstrip line coincide with 50 Q, the microwave emitted from the network analyzer 
propagates along the semi rigid cable and microstrip line, reflects at the end of microstrip 
line, and then goes back to the network analyzer. The sample is inserted between the plates 
of microstrip line as shown in Figs. 5a and b. The AC magnetic field of microwave is 
polarized perpendicular to the propagation direction. The examples of observed reflection 
coefficient spectra Su arc shown in Fig. 5c. The gradual decrease of background is caused 
by the loss of the semi rigid cable and microstrip line. The sharp dips in the 800 Oe spectrum 
are induced by the microwave absorption due to the magnetic excitations of sample. On 
the other hand, the frequency of magnetic excitations becomes higher than the measured 
frequency range at 5000 Oe. We have used the high field spectrum as the back ground and 
obtained the microwave absorption spectra ASn by the subtraction in the unit of dB. As 
exemplified in Fig. 5d, smooth spectra are obtained using this procedure even without any 
numerical smoothing. 

1. Pozar, D. M., Microwave Engineering (John Wiley & Sons, Inc. 2005). 
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FIG. 5: Microwave absorption experiment a,b, Illustrations of experimental setup and mi- 
crostrip line, c, Example of observed reflection coefficient spectra 5ii. d, Obtained microwave 
absorption spectra AS*!! by the subtraction. 
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